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SUMMARY 


Computations are presented describing the mixing and combustion of swirling 
jets in a coaxial stream. It is demonstrated that the boundary layer equa- 
tions represent the flow reasonably well until reversed flow is imminent. 

An isothermal jet with initial solid body rotation and uniform axial velocity 
yield near region profiles of axial velocity that exhibit a maximum shifted 
from the axis for swirls less than that required for the onset of reversed 
flow. For the range of parameters investigated indications are that the edge 
velocity has little effect on the behavior of the flow. Furthermore, con- 
fining the flow with a constant pressure wall, or impressing a favorable 
pressure gradient on the coaxial flow, acts to reduce the severity of the 
centerline adverse pressure gradient created by the swirl decay. 

A simple scalar eddy viscosity model, including a potential core formulation, 
is shown to describe the behavior of -weak swirling flow in the far region but 
is only in fair agreement with observations in the near region. A slight 
modification of the model, employing radial uniformity of the eddy viscosity, 
is found to be sometimes advisable. 

The effects of swirl on a burning hydrocarbon jet exhausting into a cold co- 
axial stream are shown to be intensified by the reduction of the jet density 
due to combustion. The enhanced mixing properties of high swirl flow produce 
rapid diffusion of the burning gases into the cold edge flow causing early 
cessation of the NO producing reactions. Computations show that doubling the 
initial jet swirl could reduce the NO production by 25%. 


INTRODUCTION 


There are two important features that distinguish the turbulent mixing of 
swirling streams from the mixing of non-rotating flows. The streams mix 
in a distance that can be different from that associated with the cor- 
responding non-rotating streams and may be thought of as resulting from 
an effective change in the system eddy viscosity. Furthermore, swirling 
jets with a sufficient flux of angular momentum may experience adverse 
axial pressure gradients severe enough to cause reversed (recirculating) 
flow. The scope of this study is limited to eddy viscosity effects and 
adverse pressure gradient effects only up to the point at which recircula- 
tion is first encountered. 


Because the mixing processes are so important in combustion situations it 
is not surprising that many of the observations on the effects of swirl 
have been with respect to the improvement of burning characteristics. 

Schwartz (Reference 1) demonstrated that rotating a propane-air mixture in 
an annular combustion chamber and exhausting to ambient air results in ap- 
preciable shortening of the flame length, increased flame divergence, im- 
proved stability characteristics and delayed blow-off. 

Chervinsky (Reference 2) analyzed a turbulent diffusion flame with swirl. 

He also found that the flame length was reduced and the flame stabilized 
when the fuel jet swirl was increased. Chervinsky deduced that the addi- 
tion of swirl increased the eddy viscosity and altered the turbulent 
Prandtl number. 

The set of experiments carried out by Fejer, Hermann and Torda (Reference 3) 
have provided a description of the mixing region velocity profiles for a 
swirling air jet exhausting into a non-rotating coaxial airstream under 
isothermal conditions. This data shows improved mixing due to rotation and 
is particularly helpful since zero swirl cases are presented which allow 
direct comparison between swirling and non-swirling flows. The data of 
Reference (3) also includes regimes of recirculating flow and indicates the 
existence of jet exit plane profile distortion even prior to the establish- 
ment of reversed flow. When recirculation appeared reversed flow was detected 
at the jet exit plane. Such regions are of obvious importance, in flame 
stabilization and have been studied in some detail by Syred, Chigier and Beer 
(Reference 4). 


Beer and Chigier (Reference 5) have shown that when a turbulent diffusion 
flame is confined within a rotating environment then the flame length is in- 
creased, the flame speed is reduced and the boundary of the flame is 1am- 
inarized. It was concluded that the interaction of density gradients with 
the rotating field was responsible for this behavior and for a modified 
Richardson number greater than unity the turbulence of the system can be 
damped. The definition of the modified Richardson number, R* is given in 
terms of local parameters by 


1_ 3fi. w 
p 3r r 


(— ) 
v ar' 
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This formulation assumes an invariant rotational force field. Li 1 ley 
(Reference 6) defines the Richardson number as a function of the radial 
angular momentum gradient such that a Prandtl type mixing length is in- 
creased where the gradient is negative and decreased when the gradient is 
positive. For constant density flows this behavior is consistent with the 
formulation of Rube! (Reference 7). 

In that model increased mixing due to swirl is attributed to a mechanism 
of rotational instability governed by the local criterion 

> stable 

|p (pr 4 a 2 ) = 0 neutral 

< unstable 

A correction factor for the eddy viscosity of such flows in the form of 
the pre-multiplier (1 + 90/N§ )% was used to quantitatively describe the 
relationship between mixing and swirl. 

The purpose of this report is to present some computations describing swirl- 
ing jet mixing and combustion. A jet with uniform axial velocity and solid 
body rotation exhausting into a non-rotating coaxial stream is examined as 
the jet rotation is increased until reversed flow is about to begin. The 
value of the swirl ratio, S*, at separation as well as the general nature 
of the flow in the moderate swirl regime are in agreement with observations 
(References 3 and 8). I 

In addition, the simple eddy viscosity model of Reference (7), modified to 
account for a potential core, is shown to give fair agreement with the weak 
swirl data of Reference (3) in the near region and good agreement in the 
far region (x/D > 15). Finally, finite rate chemistry calculations are per- 
formed to demonstrate a manner in which swirling flow can be utilized to 
reduce pollutants generated by the combustion process. In the case presented, 
the swirl induced rapid mixing of combustion gases with the cooler coaxial air 
flow reduces the NO production by 25%. The geometry used for these calculations 
has the basic features of the swirl -can combustor currently under development 
at NASA Lewis (Reference 9). 
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GOVERNING EQUATIONS 


The problem being described here is a turbulent swirling jet exhausting 
into a coaxial non-rotating stream. The stream can be confined or free 
and the possibility of combustion exists. The governing equations used 
are: 

Conservation of Global Mass 


fp (pur) + |p ( P vr) = 0 (1) 

Conservation of Species Mass 

p 3x p sr ' r sr W N pr r ar * w k’ k_1 ' n 
Conservation of Axial Momentum 


pu 




3 I + iL (u r 9U) 
ax r ar VP T ar ; 


(3) 


Conservation of Radial Momentum 

pr« 2 • ff (4) 

Conservation of Angular Momentum 


,, , u 3ft , 9 pvq _ 1 a / 3 an* 

P u — + P v — + 2 ^ 3 7F (y T 


ax K ar 
Conservation of Energy 


ar' 


£ p pp & + £ p 0* W ■ u If + v $ + + r2 <fF> 2) 

u . la_ r S_ r /£L rft 2 N Le 3T „ r 3a k 
k k h k r ar [ N pr y T W " c ^ U T N pr ar k C P k ar 


(5) 


( 6 ) 


The above equations result from making the boundary layer assumptions (e.g.. 
Reference 10). The invoking of this simplification allows for the solution 
of the equations by numerical methods applicable to parabolic type partial 
differential equations. At the same time the use of such methods excludes 
the possibility of calculating recirculating flows. 


These governing equations will be applied to situations that are very nearly 
recirculating. In these instances it is assumed that the radial penetration 
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of the conditions at the jet edge are primarily responsible for the behavior 
of the downstream flow field and that there is little influence from the 
downstream boundary conditions. 

The turbulent transport properties of the flow are assumed to be functions of 
a scalar eddy viscosity, yj, and the turbulent Lewis and Prandtl numbers. 

The eddy viscosity pre-multiplier model of Rubel is used in conjunction with 
the potential core formulation of Kleinstein (Reference 11), so that 


where 


W T = 0.00075 F s IPgUg-pjUj |x 
P T = 0.0286 F s r h |p e V«>(j_ u <j_l 


(core) 

(downstream) 


F s (x,r) 


F s (x,r) 


1 for |p (pr 4 a 2 ) > 0 


1, 


(oA r -‘b, 


[1 + 90 (i) for |- (prV) < 0. 

0 avg. (U.-UJ 2 3r 


J e' 


(7a) 

(7b) 


The core model is used until it matches the value of eddy viscosity given by 
the downstream model. The matching point occurs in a distance of from ten to 
twenty jet diameters and thereafter, the downstream model is used. In all 
the computations performed the Prandtl and Lewis numbers were taken to be 
unity. 

The governing equations are cast into an (x,^) coordinate system by applica- 
tion of the von Mises transformation where 


pur = ip (|£) , P vr = - ip (||) 


( 8 ) 


This transformation identically satisfies Equation (1) and recasts Equations 
(2) - (6) as 


= i 3_ , N le pur 2 3 V *k 

3X ip ^T Np^ ip dip pU 


k=l, n 


2 2 

au . L12. + lL / pur 9u \ + vrn 

ax ~ “ pu ax ip dp ^ y T \p dip ' u 2 


2 

_ ap 

u " a* 

4 

10. = . 2vn + i a , eur_ 10.) 
ax ru ~~dj W 

r c ip y 


( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 
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(13) 


r £L - IlE + pur 2 r /3Ji\ 2 . J- . IVje 3T _ r 3a k, 
C p 8X p 3x 2 y T 34> 34> Np r 34; [ C p k 34- ] 


_ 1_ z W h +IL fu IeL f£ML £L _ I£_\ i 

pu £ w k n k 4 - dip y T Np r 1 4> 3* c J 
which, along with the perfect gas law 

a k 

P = pRT z if ( 

k k 

completes the set of governing equations. 

When the boundary conditions of a. , u, a, T symmetry at the axis are applied 
to Equations (9), (10), (12) and K (13) the 1 ' Hospital rule is invoked, then, 
at the jet centerline 


da k „ N Le 3 °k , *k 
3X " “T N Pr 3 * 2 


777 * k=l, n 


3U 1 3D, , 3^U 


30 _ 0 3 (pU ) _j_ k 3 0 

'■vv ^ 1 1 v 0 


ax pu ax 


T 77 


L p 3x , )i dU ‘ w k "k 2 Np r Sf 2 
The initial and boundary conditions for the coaxial jet are 


x=0: a k -a k (4>), u=Uj(4-), fi=fij(4»), T=Tj(4»): 
J 


r or p 
w H w 

. _q. 8a k _ 3U _ 30 _ 3T _ Q 

W 34; 34> 3\J; 34- 

ip^e : a k =a k , U=u , 0=0, T=T 


P(^w } = p w (x) 


. k _ 3u _ 30 _ n 3 T _ ipn ,, , s 

* V 34< " 34 34; " * 34; * r ^w^ ~ r w^ 
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Here, the initial jet profiles and edge conditions are to be specified in- 
cluding either the pressure distribution p w (x), on the streamline originally 
located at r w or the location of the streamline, r w (x), with initial pressure 
p w . Until the wall begins to interfere with the jet inviscid relationships 
govern the edge conditions. When the wall begins to interfere with the jet 
then zero net transport across the wall is assumed. 

The governing system of partial differential equations was solved numerically 
by the explicit finite-difference technique of Zeiberg and Bleich (Reference 
12) altered to consider the presence of angular momentum and radial pressure 
gradients. The forward marching procedure is begun at the wall or edge 
streamline where the axial pressure gradient is known (prescribed or assumed). 
In order to continue this scheme the axial pressure gradient along the neigh- 
boring streamline must be obtained. This is found by inward integration of 
the radial momentum equation. The technique can now be continued along each 
streamline to the axis. Where ,a value for the radial velocity is required it 
is computed from the streamline deflection, i.e., v = u(3r/3xL, using back- 
ward differences. At the start of the calculation the radial velocity is 
assumed to be zero. 

In the case where the edge streamline position is prescribed the associated 
pressure must be iterated upon for each axial step. The iteration is com- 
pleted when the assumed pressure results in a streamline position that is in 
agreement with the prescribed position to within one percent of the radial 
mesh width. 

The basic concepts of Ferri , Moretti and Slutsky (Reference 13) have been ap- 
plied to numerically couple the mixing and finite rate chemistry equations. 


7 



RESULTS AND DISCUSSION 


A computer program developed at ATL based upon the analysis given above 
was used to determine the effect of several parameters on the behavior of 
a swirling jet exhausting into a coaxial stream and to explore application 
of these principles to combustor studies. 

Isothermal Jet With Uniform Initial Conditions - This example was 
chosen so that the basic behavior of a swirling jet could be studied and 
the influence of various parameters ascertained. A 2.54 cm diameter jet 
with uniform axial velocity of 122 m/sec and solid body rotation is assumed 
to be exhausting into a 31 m/sec non-rotating coaxial stream bounded by a 
15.24 cm diameter wall. The pressure and temperature of the streams are 
1 atmosphere and 300K respectively. These conditions have integral properties 
similar to those of Reference (3). 

The effect of swirl on the axial velocity decay at the centerline is shown in 
Figure (1). An adverse pressure gradient, created by the decaying swirl field, 
slows the flow increasingly as a function of initial swirl up to Rflj/Uj = 49/48. 
Further increases in swirl caused failure of the program due to reversed flow. 
Increased swirl enhances mixing since the eddy viscosity is greater (unstable 
angular momentum distribution) and,, for Up/Uj < 1, the reduced centerline axial 
velocity is also in the direction of the final fully mixed solution. 

In Figure (2) the streamlines are plotted for the Rfij/Uj = 49/48 case (note 
stretched radial coordinate) and demonstrate the deflection of streamlines 
near the region of reduced axial velocity. One can anticipate the separation 
zone that will appear upon further increase of the deflection of the stream- 
lines near the axis. 

For the initial profiles and eddy viscosity model selected here the effects 
of adverse pressure gradient crest at x/D % 5. The radial distribution of 
axial velocity at this location is shown in Figure (3) as a function of 
Rftj/Uj. The effects of increased rates of mixing extend the mixing zone to 
greater radial distances with increased swirl. At the same time the adverse 
axial pressure gradient causes the maximum velocity to be shifted away from 
the axis for Rq./U, > 3/4. 

J J 

The downstream development of the axial velocity distribution is given in 
Figure (4) for Rfij/Uj = 1. The maximum velocity has been shifted from the 
centerline within J one diameter of the jet exit plane. The bulge in the pro- 
file grows until x/D % 5 where the shear forces begin to dominate the inner 
region and eventually return the profile to that of ordinary mixing (x/D £ 10). 

The distribution of angular velocity with downstream distance for Rfij/U,- = 1 
(Figure 5) also has a peak away from the axis. Shear acts to reduce J the 
angular velocity near the interface of the jet and edge stream but, in the 
nearly inviscid central region of the flow (x/D > 5) the angular velocity is 
also reduced. In this region the reduced axial velocity has caused the 
streamlines to be deflected outward (Figure 2). Since angular momentum (i.e.. 
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U # ■ 31 m/sec 



FIGURE 1. EFFECT OF SWIRL ON CENTERLINE AXIAL VELOCITY DECAY 
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FIGURE 5. DECAY OF RADIAL DISTRIBUTION OF 
ANGULAR VELOCITY 
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p 

r n) along a streamline remains constant under these circumstances, then the 
angular velocity, Q, along streamlines near the axis must be decreasing. The 
effects of shear begin to dominate the swirl behavior near the jet axis at 
x/D % 5 as was the case with the axial velocity. 

The difference in pressure between the jet edge and centerline is another 
indication of the concentration of angular momentum within the jet. Figure 
(6) gives the axial decay of this pressure difference as a function of ini- 
tial swirl. It appears that within five diameters of the jet exit plane there 
is little remaining of the initial radial pressure gradient. Since the length 
for decay is not much dependent upon the initial swirl, it follows that for 
nearly recirculating flows, the adverse pressure gradient impressed upon the 
flow at the axis is directly proportional to the initial pressure difference 
between the jet edge and centerline. 

From these calculations the near flow field of a swirling jet exhausting into 
a coaxial stream can be described. Initially there is mixing at the inter- 
face of the jet and edge streams that acts to diffuse angular momentum from 
the jet. This diffusion acts to increase the pressure and, thus, to reduce 
the axial velocity throughout the jet. The streamlines move away from the 
centerline with the greatest deflections near the axis. This acts to further 
reduce the angular velocity near the jet centerline causing the most severe 
adverse pressure gradients to be located along the axis. It is for this 
reason that the axial velocity profile has its maximum shifted away from the 
axis. If the flow is not recirculating then shear gradients will begin to 
act at the centerline to increase the velocity before a stagnation point is 
attained, eventually producing the usual downstream mixing profile. 

The effect of an exterior favorable pressure gradient on the behavior of the 
jet is shown in Figure (7). As expected, such a free stream gradient acts to 
reduce the centerline adverse pressure gradient. In the case presented here 
(i.e., p w = 1-0.01575 x/D atmosphere) little velocity decrease is noted at 
the axis for initial swirl Rftj/Uj = 1 and no reversed flow was encountered 
even for swirl as great as Rfti/Uj = 5/4. The peak effects of swirl still 
appear to be restricted to within the first five jet diameters. 

The effects of confining the jet are shown in Figure (8). Here the pressure 
at the confining wall was maintained constant and the initial wall radius 
was only twice the jet radius. The velocity profiles are unchanged (see 
Figures 4 and 5) until the wall interferes with the flow at x/D > 2. The 
wall reduces the diffusion of angular velocity which acts to delay the decay 
of the radial pressure gradient. The axial pressure gradient is made less 
severe so that the minimum velocity at the axis is greater than for the un- 
confined flow. Of course, the downstream results exhibit higher velocity 
since the amount of coaxial flow available for mixing is limited by the con- 
finement. 

The effects of the edge velocity on the swirling jet are indicated by Figure 
(9). Increasing the edge velocity from IL/lk- = 1/4 to U e /Uj = 3/2 had little 
effect on the value of the minimum velocity at the jet axis, but it did act 
to move this point toward the jet exit plane. The adverse axial pressure 
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Reduced centerline pressure* 



FIGURE 6. EFFECT OF SWIRL ON THE DECAY OF RADIAL 
PRESSURE DIFFERENCE 
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Reduced axial centerline velocity. 



FIGURE 7. EFFECT OF SWIRL ON CENTERLINE AXIAL VELOCITY 
DECAY IN PRESENCE OF FAVORABLE EDGE PRESSURE 
GRADIENT 
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FIGURE 8. VELOCITY PROFILE DECAY WITH WALL INTERFERENCE 
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FIGURE 9. EFFECT OF EDGE VELOCITY ON CENTERLINE AXIAL VELOCITY DECAY 





gradient at the axis is increased with increased edge velocity but axis 
shear gradients balance this effect. 

Uniform initial profiles were used here as the simplest means of generating 
results for the effects of basic parameters on swirling jet behavior. These 
profiles preclude comparison with current data except, perhaps, on some 
integral basis. The characteristic parameter for describing swirling jets 
exhausting into a quiescent medium is the swirl parameter, 

f" puAdr 

5 = 2 ° 

R[ u 2 rdr - (P -P) rdr] 

J 0 Jo 


( 16 ) 


which is an invariant dimensionless ratio of angular momentum to axial momentum. 
Since S is no longer invariant when a coaxial stream is involved, it is neces- 
sary to define the parameter, S*, in the same manner as Equation (16) but with 
the jet radius, R as the outer limit of the integrals. This parameter reduces to 
S for the quiescent case but does not account for the presence of a coaxial 
stream. Since the edge velocity does not seem to influence the behavior of 
the near flow, except for slight displacement in the axial direction, this 
parameter may prove useful for prediction purposes (e.g., onset of flow reversal). 

For uniform jet profiles of axial and angular velocity and constant density it 
is easy to demonstrate that 



The numerical computations have indicated that the maximum axial velocity be- 
gins to be displaced from the axis at S* % 0.43 (Figure 3) and reversed flow 
is imminent at S* £ 0.69 (Figure 1). Chigier and Chervinsky (Reference 8) 
found these points to occur at S ^ 0.5 and l £ 0.6, respectively for a tangen- 
tial injection swirl generator. Of course, the jet exit plane profiles are 
also important in determining the downstream behavior of the flow. If the 
jet exit plane flow exhibits reduced velocity near the axis then one would 
expect flow reversal to occur at lower values of S* (calculated at the jet 
exit plane) . 

Decay Of A Swirling Jet With Nonuniform Initial Conditions - Computa- 
tions of swirling jet flow were compared with some observations of Fejer et.al., 
(Reference 3) using the measured velocity profiles as the initial conditions 
for the calculations. The centerline decays of axial velocity for S* = 0, 0.14 
(Cases la, lb of Reference 3) are presented in Figure (10) and it is seen that 
the behavior of the non-rotating case is adequately described by the eddy 
viscosity formulation of Equations (7a) and (7b). For S* = 0 the velocity pro- 
files are also well described by this model (Figure 11). 
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FIGURE 12. VELOCITY PROFILE DECAY FOR u-m(x,r), S # *0.14 


For the case of weak swirl (S* = 0.14) only a fair description of the axial 
velocity decay in the first ten diameters is given by the pre-multiplier 
eddy viscosity model but downstream of x/D = 15 the agreement is excellent. 
The velocity profiles (Figure 12) are also in agreement beyond x/D = 15 and 
compare reasonably well in the outer mixing region, even for x/D < 15, but 
in the inner region the discrepancy is apparent. 

The eddy viscosity model based on local rotational stability was modified so 
that, at any axial position x. 
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any r. 


This model implies that at a given downstream location, if any segment of the 
flow is rotationally unstable then the eddy viscosity is increased over the 
entire jet width. Calculations using this model showed somewhat improved 
axial velocity profiles but increased the discrepancies between the swirl 
velocity profiles (Figure 13). Nevertheless, the more rapid decay of swirl 
given by the numerical computations is consistent with results presented in > 
References (2) and (6). In order to fully describe both axial and swirl 
velocities a tensor eddy viscosity is inferred (e.g.. Reference 6). 

The radial behavior of eddy viscosity is especially important in nearly re- 
circulating regions since shear gradients are, in general, responsible for in- 
creasing the axial velocity near the axis. Numerical calculations with ini- 
tial conditions provided by the post-recirculation profiles of Reference (3)'s 
Case Id, have indicated some superiority of the model y=y(x) in this regime. 
Such a formulation also acts to delay the onset of reversed flow and can be 
used to bypass such regimes for some high swirl computations where downstream 
and outer region mixing are of primary importance. For this reason the pre- 
multiplier of Equation (18) is used for all further results presented here. 

Applications To Finite Rate Combustion - In order to demonstrate 
some effects of jet swirl on the finite rate chemistry process, computations 
were made for a stoichiometric hydrocarbon-air jet exhausting into a channel 
of cold air. For comparison, two sets of calculations were performed dif- 
fering only in swirl parameter (i.e.. Case 1, S* = 0.36; Case 2, S* = 0.165). 
The geometry of the problem, as well as the initial conditions, are specified 
in Figures (14) and (15) and have basic parameters representative of the swirl- 
can design of Reference (9). The jet mixture is made ignitable by assuming 
that its exit plane temperature is 1130K. The chemical reactions used for 
the computations are given in Table I and their reaction rates can be found 
in References (14), (15) and (16). 
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FIGURE 13. VELOCITY PROFILE DECAY FOR u»w(x), S # =0.14 
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FIGURE M. SCHEMATIC OF COAXIAL SWIRLING JET MIXING AND COMBUSTION 
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FIGURE 15. INITIAL PROFILES POR CASE 1 AND CASE 2 
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CHEMICAL REACTION SYSTEM 
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The temperature and velocity profiles at three axial stations for Case 1 are 
presented in Figures (16a), (16b) and (16c). At the low velocities used for 
these calculations combustion of the jet mixture is very rapid and at x = 3 cm 
the axis temperature has risen to 2700K. The swirl and axial velocities have 
profiles that are similar to those previously presented and show typical axial 
decays toward uniformity. The axial velocity profile exhibits a maximum away 
from the centerline even though S* < 0.43. This is a result of the reduced 
density in the burning region. The adverse pressure gradient established 
by the swirl decay, when acting on fluid with reduced inertia, tends to pro- 
duce the same flow patterns as in the isothermal cases, but at reduced values 
of S*. 

The mass fractions of the more important species are presented in Figures (17a), 
(17b) and (17c). It is to be noted that NO is being produced through much of 
the jet region, where the temperature is high. The local maximum of methane 
mass fraction appearing near the edge of the mixing region at 3 cm is due 
to methane that has diffused to a region of lower temperature wnere its combus- 
tion has been delayed. 
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FIGURE 16. CASE 1 VELOCITY AND TEMPERATURE PROFILES (S*-0.36) 
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The temperature and velocity profiles for Case 2, are presented in Figures 
(18a), (18b) and (18c). The swirl profiles are similar to those of Case 1 
but reduced in scale. An axial velocity maximum of small proportion occurs 
away from the axis although the value of the swirl ratio is only S* = 0.165. 

The species distributions for Case 2 (Figures 19a, 19b and 19c) show behavior 
little different from those of Case 1 and are included for completeness. 

The effect of swirl on overall pollutant production is shown in Figure (20). 

Here the emissions of CO and NO are shown for Cases 1 and 2 as a function of 
axial distance. The CO production curves are almost identical for the two 
cases. Burning takes place earlier for the high swirl case because of the 
reduced jet velocity and, therefore, the CO production curve for S* = 0.165 
slightly lags that for S* = 0.36. The behavior of the NO production curves is 
somewhat different. Until x £ 4 cm there is more NO produced by the high 
swirl case due to the longer residence times. The reactions producing NO, 
however, are relatively slow and highly temperature dependent. Thus, the 
increased mixing due to increased swirl more rapidly diffuses the burned gases 
into the cold edge stream causing an earlier cessation of NO production. For 
the two cases calculated here the NO production was reduced by 25% upon doubling 
the initial swirl. 
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FIGURE 18. CASE 2 VELOCITY AND TEMPERATURE PROFILES 
($*=0,165) 
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FIGURE 19. CASE 2 CHEMICAL SPECIES PROFILES (S*=0.165) 
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FIGURE 20. EFFECT OF SWIRL ON POLLUTANT PRODUCTION 
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CONCLUSIONS 


Computations have been presented describing the mixing and combustion of 
swirling jets in a coaxial stream. It is demonstrated that the boundary 
layer equations represent the flow reasonably well until reversed flow is 
imminent. 

For an isothermal jet with initial solid body rotation and uniform axial 
velocity the profiles of axial velocity exhibit a maximum that is shifted 
from the axis for swirl ratios S* > 0.43. The onset of reversed flow is 
found at S* % 0.69. For 0.25 < U e /Uj < 1.5 the edge velocity has little 
effect on the behavior of the flow. J The calculations indicate that con- 
fining the flow with a constant pressure wall, or impressing a favorable 
pressure gradient on the coaxial flow, acts to reduce the severity of the 
centerline adverse pressure gradient created by the swirl decay. 

The simple scalar eddy viscosity model of Rubel , including a potential core 
formulation, is shown to describe the behavior of weak swirling flow in the 
region x/D > 15 but is only in fair agreement with the data of Fejer et.al. 
in the near region. A slight modification of the model, employing radial 
uniformity of the eddy viscosity, is found to be sometimes advisable. 

The effects of swirl on a burning hydrocarbon jet exhausting into a cold 
coaxial stream are shown to be intensified by the reduction of the jet 
density due to combustion. The enhanced mixing properties of high swirl 
flow produce rapid diffusion of the burning gases into the cold edge flow 
causing early cessation of the NO producing reactions. Computations show 
that doubling the initial jet swirl could reduce the NO production by 25%. 
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